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Abstract 

In this note we present a new numerical method for solving Lattice Quantum Field Theory. 
This Source Galerkin Method is fundamentally different in concept and application from Monte 
Carlo based methods which have been the primary mode of numerical solution in Quantum Field 
Theory. Source Galerkin is not probabilistic and treats fermions and bosons in an equivalent 
manner. 
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1 Introduction 



Over a decade ago, it was realized that the algorithms and computer power existed to calculate non- 
perturbative consequences of Quantum Field Theory 0] . Initial results using Monte Carlo methods 
and the quenched approximation were a striking confirmation of the power of numerical methods 
and the general correctness of the ideas of lattice QCD ||. After countless hours of calculation 
with ever more powerful computers and vast improvements of the initial class of algorithms, there is 
still the hope of great success as machines become much faster. While there has been demonstrable 
progress, it can be argued that important but mostly incremental improvements in results have been 
made. 

It is the purpose of this note to introduce a very different class of approaches to numerical 
quantum field theory. The new methods, while not without their own difficulties, have the promise of 
allowing a great increase in the accuracy and speed of numerical calculation. Particularly important 
is that fermions can be treated in essentially the same manner as bosons. 

For clarity, we will illustrate these ideas with very simple examples. Our initial analysis was 
produced through the use of ALJABR* (a product of Fort Pond Research) which is a very reliable 
computer algebra program. It is not necessary to resort to Fortran or C coded implementations 
of our methods until quite complicated systems are studied. A good workstation will produce a 
remarkable number of results. 

Our approach is based on considering field theory in the presence of external sources. We start 
with the functional equations for the vacuum amplitude Z . On a lattice, these are converted to a 
set of coupled linear differential equations in the discretized sources J. 

The high degree of symmetry due to the translation and reflection invariance of the lattice can 
be used to our advantage in constructing the solution. For this note, Z will be approximated by a 
power series in J |3| . Faster convergence and the full power of our ideas is obtained with other more 
structured approximations that will be presented elsewhere. 

Truncated approximations in multiple dimensions are intrinsically inconsistent. One particularly 
elegant method to handle this problem is to require that weighted averages of the equation residual 
vanish. A clever choice of weight functions speeds the convergence to the exact solution. Techniques 
of this sort are familiar in the numerical solution of differential equations and are called Galerkin 
methods The accuracy of the approximation can be judged by measuring its numerical stability 
and convergence. 



1 



Since our expansion and weight functions depend on the external sources, we call this method 
the Source Galerkin method. Except for the anti-commutation of fermionic sources, fermions and 
bosons can be treated symmetrically. Nested solutions for dynamical fermion problems become 
possible with existing computers. 



2 Functional Formulation in the Presence of Sources 

We begin by studying a scalar field $(ar) in the presence of an external source J(x). As an exam- 
ple, consider a quartically self interacting scalar field $ on a periodic D-dimensional lattice, with 
Euclidean action 

s E (^j)=J2lT,{^(.- a + M2 ) x ,y^y)+^ 4 ^ J ^ > C 1 ) 

x y 

where 

(-□ + M 2 ) x y ee (2D + M 2 ) S x>y - <W • (2) 

A 

The lattice spacing a has been absorbed in the definition of the sources and parameters, and the p, 
summation is over nearest neighbors. The generating function Z(J) satisfies the the following set of 
coupled partial differential equations in source space 

There is one equation for each J x . The Euclidean (unconnected) lattice Green's functions can be 
extracted from Z(J) by differentiation 

d m Z 

$(x m )) = G m (xi,... x m ) = — k-— to = 1,2,3... . (4) 

O J Xl . . . OJ Xm J=0 

For the rest of this note we will examine only D = 1 theories; the ideas generalize to higher 
dimensions. If fermions were present, the Grassman nature of the equations causes only technical 
complications. 

3 Boundary Conditions 

A straightforward way to attack (||) is to construct a power series solution |jj 
Z(Ji, J%, ■ ■ • Jn) = G ni . n2t ... nN 

J? J?... J™ , (5) 



ni,n 2 ,...,njv 



with Hi, i = 1,2, . . . ,N nonnegative integers. 

Before solving (||) for the coefficients G nitn2t ,,, nN , we must specify boundary conditions. Since (||) 
is homogeneous in Z, we can normalize the vacuum amplitude 

Z{J = 0) = 1 , (6) 

and require odd Green's functions to vanish 

3 n 7 

= in WT~ =0 VJ„n= 1,3,5... . (7) 

d J Kl . . . 9 j=o 

If equation (^|) is considered as an initial value problem, we need also to specify the second derivatives 
at the origin 

87 

This appears to require a prior knowledge of the two point Green's function j^] . 

Alternatively, we demand the constructed 7 to approach smoothly the free field generating 
functional 

Z (J) = e J ^- D + M ^.y J y , (9) 

when g — ► + . This requirement is an implicit boundary condition of all Monte-Carlo (path inte- 
gral) simulations and, together with the above conditions, uniquely determines Z. Other classes of 
solutions to the field equations (||, JlO| , J?]]) may also be of physical interest. 

We found that one way of implementing this boundary condition is to truncate Z at some finite 
order K in the sources J, so that all unconnected Green's functions with more than K legs vanish. 
Truncation then introduces an approximation scheme which can be systematically improved by 
truncating at successively higher orders and taking the limit K — > oo. The convergence of this 
process will be discussed in future publications. 



4 Lattice Symmetries 

The number of coefficients in (j^) grows exponentially with the number of sites and with increasing 
truncation order j^] . Consequently, this particular version of the Source Galerkin method is only of 
use for small lattices. We can get somewhat further by exploiting symmetries to reduce the number 
of unknowns to solve. For instance, in a lattice of ten sites, one can invoke translation and reflection 
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symmetries and group the 45 second order monomials of the form J X J V into five invariant classes, 
thus reducing considerably the number of unknowns (see (|i9|)). 

Since the differential operators in ([}]) are closed under the action of the lattice symmetry group, 
it is natural to look only for invariant solutions Z . 

a n 

where Pi n \j) are invariant polynomials of order n corresponding to a particular symmetry class a. 
The number of invariant classes for a given order will depend on both the number of lattice sites 
and on the number of symmetry operations for a given lattice. Details of how to construct these 
invariant polynomials will be given elsewhere, but it is important to realize now that they form a 
complete set for the class of lattice invariant solutions we are considering. We believe the ability to 
easily exploit symmetries is a major advantage of our approach. 

Once Z is constructed as lattice symmetric, it is sufficient to use only one equation from the 
coupled set ([5j): all equations in the coupled system express identical dynamics and only differ by 
lattice translations and reflections. 



5 The Galerkin method 

When solving a differential equation by power series methods, we substitute the series into the 
equation. Equating like terms in the resulting polynomial yields a system of linear algebraic equations 
that determine the expansion coefficients of the original power series. 

Applying this procedure to (^) and the power series (0) for Z yields an inconsistent set of linear 
equations for the truncated Green's functions: there are more independent equations than unknowns 

. The Galerkin method is particularly suited to deal with this problem. 

The key features of the Galerkin method can be stated concisely [Q. Assume that we have a 
problem governed by a linear differential equation 

Lf = b . (fl) 

L is a linear differential operator, and b a known function of the N variables x = (xi, . . . ,Xn). 
The solution is sought in a N-dimensional domain T>. The Galerkin method assumes that the exact 
solution / = /(x) is accurately represented by an approximate 

K 

/A'( a ; x ) = <Pa(x) +^aj<Pj(x) , (12) 

3=1 
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where tpo is introduced to satisfy the boundary conditions and the aj's are coefficients to be deter- 
mined. The K functions tp are selected according to some criterion, like knowledge of the approximate 
behavior of the solution or mathematical completeness. These are called basis, or trial, functions. 
Substitution of (|lj) in ( |TT| ) produces a nonzero residual R given by 

K 

i?(a;x) = L(f) = LM+^ajL&j) . (13) 

3=1 

The unknown coefficients aj are determined by solving the following system of linear algebraic 
equations 



< tp k ,R > = < tp k ,b > , (14) 



K 



— ► < PkiLcpj > aj —< fk,b > - < (fik,LLp > , k = l,2, ...K , 

3=1 

where Cp are suitable test functions, and < , > is the inner (scalar) product 

< 0,X >= [ d N x<j> X ■ (15) 
Jv 

In the pure Galerkin method, the test functions are chosen from the same functional space as the 
basis functions <p in (|l^). 

Under very general conditions Q| it is possible to prove that the Galerkin solution converges 
to the true solution / in the mean 

Jim / H/;,-- /|| =0 . (16) 



6 Example 

As an illustration of the Galerkin method with truncated boundary conditions, we solve ([3|) for 
the case of a one dimensional lattice with three sites (D = 1 and N = 3 in (||)), expanding the 
generating function Z up to the first non-trivial order in the sources, that is, up to quartic invariant 
polynomials. 

First, we construct Z\ = Zl(Ji, J2, J3) as a linear combination of lattice invariant polynomials. 

zkji, j 2 , j 3 ) = 1 + 4 2) p (2) (j) + a ^p[ 2 \j) + 4 4) p (4) (j) • ■ ■ + 4 4) ^ ( V) . (17) 

with 

p^(j) = j 2 1 +j! + j! , (is) 



P X (2) (J) = JlJ 2 + J 2 J 3 + J 3 Jl 

pf\J) = Jt + Ji+J$ 



(n) 

Only six P„ are needed. The boundary conditions in section 3 have already been imposed. We 
choose as our integration domain T> = [— e, e] 3 , a hypercube of size 2e centered at the origin, with e 
infinitesimal and set to zero after the calculation is done (this is of course possible only with algebraic 
languages like ALJABR, otherwise a numerical limit has to be taken). 

Because the operator in (|^) is odd and T> is symmetric around J = 0, we need odd test functions 
in order to obtain non-vanishing inner products. A prescription for test functions that we found 
numerically convenient, since it provides the necessary number of independent equations for the 
coefficients ai™ , is the derivative of the basis polynomials respect to a source J x 

apt") 

r<"> = ^- . (19) 

Lattice symmetry makes any choice of x equivalent. We set x = 1. Any other set of test functions 
that produces the correct number of independent equations will also work, but with a different rate 
of convergence. 

The Galerkin method is then applied to the residual of (||). The resulting coefficients aS 
are rational functions of the mass and the coupling. For instance, the Galerkin approximation to 
($(0)$(0)) is 



<*(0)«(0)>- = =2a?) S , (20) 



dZ% 

dJf j=o~^° ~ B(M, g) 



with 



A(M, g) = M 2 (M 2 + 1)(M 2 + 3)(266M 4 + 929M 2 + 788) + (21) 
35g(M 2 + l)(M 2 + 2)(14Af 2 + 27) , 

B(M, g) = M 4 (M 2 + 3) 2 (266M 4 + 929Af 2 + 788) + (22) 
5 (1288Af 8 + 7778M 6 + 17997M 4 + 18759M 2 + 7092) + 
105 5 2 (M 2 + l)(14M 2 + 27) , 

One can check that the above expression (ptj|) has the right free field limit when g — > + . 
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The convergence of this process is illustrated in Table 1, where we show our results for the two 
point function 



G, 



d 2 Z 



dJ;dJj 



on an 11 site lattice with M = 1 and g = 0.5. We see that the results agree extraordinarily well 
with the Monte-Carlo simulation, and that the convergence is oscillatory M. 



7 Conclusions 

We have presented a new method to numerically attack problems in lattice QFT. It is not probabilis- 
tic and possesses extreme flexibility, allowing for maximum input of information (symmetry, spectral 
properties, etc.). Systematic improvements of the approximations are possible, forcing convergence 
(in the mean) to the exact solution. We have illustrated the method by looking at truncated Tay- 
lor series expansions. These lose tractability for big systems. We will present more useful nested 
approximations elsewhere. 

The differential formulation requires treatment of boundary conditions. Truncation (among 
other boundary conditions Q) guarantees that the solution obtained from the differential equations 
corresponds to the solution of the path integral. Other types of solutions can also be analyzed with 
our method. 

Finally, the Source Galerkin method gives real hope that lattice fermion theories will yield to nu- 
merical treatment. Monte Carlo approaches in QCD involving fermions are problematic. Calculation 
of the fermion determinant makes inclusion of dynamical quarks in lattice QCD prohibitive. Simula- 
tions of strongly correlated systems in condensed matter physics are afflicted by the minus sign prob- 
lem (||). We expect Source Galerkin to be a powerful tool for studying these systems since the above 
problems do not affect our method. The difference in computational complexity between bosons and 
fermions is minimal. Indeed, part of the strength of our method is the symmetry of approximations 
for bosons and fermions. 
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Tabic 1 



i-j 


J 4 


J 6 


J 8 


Shanks 


M-C 





0.326538 


0.360178 


0.347223 


0.35082 


0.3510 ± 0.0005 


1 


0.094508 


0.117673 


0.108121 


0.11091 


0.1110 ± 0.0004 


2 


0.026768 


0.038823 


0.033422 


0.03509 


0.0351 ± 0.0003 


3 


0.007454 


0.012940 


0.010269 


0.01114 


0.0112 ± 0.0004 


4 


0.002075 


0.004461 


0.003197 


0.00363 


0.0036 ± 0.0003 


5 


0.000694 


0.001922 


0.001218 


0.00147 


0.0011 ± 0.0004 



Tabic 1: Results for the two point function dj in a N=ll lattice (anharmonic oscillator). Columns: 
\i — j\, 4th., 6th. and 8th. order Galerkin approximations, numerical extrapolation (Shanks trans- 
formation) and Monte-Carlo data with statistical error. 
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